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Abstract 

The galileon model was recently proposed to locally describe a class of modified gravity theories, 
including the braneworld DGP model. We discuss spontaneous symmetry breaking of the self- 
accelerating branch in a multi-galileon theory with internal global symmetries. We show a modified 
version of Goldstone's theorem is applicable to the symmetry breaking pattern and discuss its 
implications. We also derive the Hamiltonian of a general multi-galileon theory and discuss its 
implications. 
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I. INTRODUCTION 



The DGP model [TJ [2] is a 5 dimensional braneworld theory that non-trivially modifies 
General Relativity (GR) in the infrared. Nevertheless, at sub-crossover (sub-Hubble Length) 
scales many of its properties can be captured by a 4D (boundary) effective theory [21 H]. 
This effective theory amounts to GR coupled to a scalar field ix whose equation of motion 
has only second derivatives and is invariant under the Galilean shift tt — » ix + a M x M + b, 
and b being constant. This scalar is related to the bending of the DGP brane in the bulk 
and has been termed as galileon [5]. 

As ghost instability has been identified on the phenomenologically interesting self- 
accelerating branch of the DGP model [6] , which can also be easily seen in the local galileon 
approximation [31 H] , attempts have been made to generalize the DGP galileon description to 
produce a healthy modified gravity theory [5ll7HTo]. In [5], the authors wrote down the most 
general single galileon Lagrangian. Remarkably, there are only d + 1 possible galileon terms 
in d dimensional spacetime, and ghost free self-accelerating background solutions have been 
shown to exist in a generalized galileon theory. However, a few phenomenologically challeng- 
ing problems have also been identified in the single galileon theory, such as Cherenkov-like 
radiation in the solar system, superluminal propagation far away from a matter sauce and 
very low strong coupling scales [5]. It turns out that these problems can be avoided by 
adding another galileon (in a bi-galileon theory), meaning the theory space of the single 
galileon model is actually too small [13]. A local bi-galileon description is also what one 
might expect from co-dimension 2 braneworld models [121 H3 E] , as there are generally two 
brane bending directions. 

One would want to generalize the galileon description to have even more degrees of 
freedom [TT1 [131116] . To avoid a proliferation of possible terms in the theory, we can impose 
internal (global) symmetries within the multiple galileons, so that the multiple galileons form 
some representation of a group [14J, ix = (jn, ix^). That is, the multi-galileon Lagrangian 
is imposed to be invariant under the internal transformation 

7Ti -+ TV'tT,-, (1) 

where TZ/ is the representation matrix of a certain group and summation over repeated group 
indices is implied. Notice that the internal symmetry could originate from braneworld sce- 
narios, as has been identified for the SO(N) fundamental representation [Tl"l [To]. For other 
interesting field theoretical and cosmological implications of the galileon theory, see [T8H22] . 

In [Hj, we wrote down all possible multi-galileon terms that are consistent with the 
fundamental and adjoint representations of SO(N) and SU(N), and looked for soliton solu- 
tions in multi-galileon theories; We did not consider coupling the symmetric multi-galileon 
to gravity. In this paper, we will put the symmetric multi-galileon in the context of modified 
gravity. In Section [II], we will venture a tentative coupling, but we want to emphasize that 
the main results of this paper are insensitive to this explicit coupling. In Section III we 
discuss the spontaneous symmetry breaking phenomenon of the symmetric multi-galileon 
theory on a self-accelerating background. Starting from an example, we build up a new 
version of Goldstone's theorem in symmetric multi-galileon theories that for every broken 



continuous symmetry a canonical kinetic degree of freedom is lost. In Section \P7\ we derive 
the Hamiltonian formulation of a general multi-galileon theory (with or without internal 
symmetry) and find it is not bounded below. We speculate whether this might be overcome 
in more complete theories. 
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II. MULTI-GALILEON MODIFIED GRAVITY 



In the original galileon model [5], the galileon is coupled to graviton mainly via the kinetic 
mixing 

hfj, v = + 2tt^, (2) 

where h^ v and are Jordan and Einstein frame (perturbative) metrics; 7r's contribution to 
the energy momentum tensor, or, its direct influence to the geometry is negligible. So in a 
sense the galileon modified gravity is a "genuine" infrared modification of General Relativity, 
differing from models such as quintessence [23], which has significant contribution to the 
energy momentum tensor. In this paper, we stick to this paradigm and tentatively propose 
the multi-galileon's coupling to gravity as 



S = d 4 x 



M 2 1 ~ 

■-^/i^V + 2^ + fa + - + **) T + C « 



(3) 



where T = i] fJil/ T flu and C n is the multi-galileon Lagrangian. For a general multi-galileon 
theory without internal symmetries, we might want to redefine 7t[ = ~K\ + ... + to simplify 
the coupling, while keep the structure of C w unchanged. But this is usually not feasible in 
symmetric multi-galileon models. For example, in the case of SO(N) fundamental represen- 
tation, 7r = (tti, 7T2, 7Tjv) can not be linked to n' = [tt[, 7T2, ttn) by an internal SO(N) 
transformation. (Note that the SO(N) invariant coupling P(n 2 )T, P{^ 2 ) being a general 
function of 7rVj, has been considered in [16], and the authors found gradient instability as 
well as superluminal excitations for the spherically symmetric background.) We could argue 
that from the viewpoint of braneworld scenarios the coupling (|3]) (instead of, say, itiT) might 
be what one might expect for symmetric multi-galileon models. In a braneworld setup, the 
multiple galileon fields living on a brane usually descend from the extra dimension coordi- 
nates as functions of the 4D brane volume coordinates [HI US]- Since the symmetric multiple 
galileon fields enjoy some internal symmetry, the extra dimensional coordinates must have 
the corresponding symmetry at least near the brane. As the near brane geometry is expected 
to plays a role in determining the coupling to gravity, we may expect the different multiple 
galileons couple to gravity on a equal or similar basis. 

At distances and time scales shorter than the Hubble length, the Friedmann-Robertson- 
Walker metric can be considered as a perturbation above Minkowski spacetime. Due to the 
kinetic mixing @, the cosmic profile of the multi-galileon can be cast within the Hubble 
length as [5] 

Stt = -h,H 2 - Hl)x^ + 1 -{H - H gI )t 2 , (4) 

where S7r = 7Ti + ... + tc^, H is the actual Hubble parameter for a given source and 
H gr is the hypothetical Hubble parameter in GR with the same as the source. Thus the 
cosmic background configuration of S7r is given by — \{H 2 — H^x^x^. Assuming all the 
fields have similar coordinate dependence, the vacuum solution is given by 

7f< = -^hx^, T l k = k 1 + ... + k N = H 2 -H 2 I . (5) 



3 



III. GOLDSTONE'S THEOREM IN SYMMETRIC MULTI-GALILEON MODI- 
FIED GRAVITY 

In this section we will see that the symmetric multi-galileon modified gravity exhibits 
spontaneous breaking of symmetries on some vacuum solution, and for every broken con- 
tinuous symmetry the theory loses a canonical kinetic term, which resembles the usual 
Goldstone's theorem in a scalar field theory. We will also discuss the implications of this 
modified Goldstone's theorem. 



A. An Example 

Let us first see a simple example of this theorem: spontaneous breaking of the SO(N) 
fundamental representation. The most general SO(N) multi-galileon Lagrangian in the 
fundamental representation is given by [H] 

= -a d^d^ - (3 6^ ] d p Tr i d% i d^d x d\ h (6) 

where 8yJ^ = 3\6? a 6%6^, and a and (3 are free parameters. Varying (3) with respect to 7Tj, 
we get the equations of motion: 

2a Dtt, + 4/3 6^,8^3^0x0^ = -T. (7) 

We would like to see whether there is any self-accelerating background (or vacuum) in this 
theory. By a self-accelerating background, we refer to the case where the universe has a (at 
least approximately) de Sitter solution without support of a cosmological constant, i.e., the 
case where 7fj = — ^kiX^ with E/c = H 2 > and T = is a solution to the equations of 
motion ([7]). Substituting this profile into the equations of motion, we get 

- Ak(a + 30kj) = 0, (8) 

which reduce to 

h = 0, (9) 

or 

vh = do) 

The former solution corresponds to Minkowski spacetime, while the later can be a self- 
accelerating solution if a//3 < and Hk = H 2 > 0, which we assume to be satisfied. Note 



that (10) is not an isolated solution, instead it is a continuum of possible solutions. 

Then we would like to see whether the self-accelerating solution can be free of ghosts, 
negative canonical kinetic terms. To this end, we expand the Lagrangian (|6| above the 



background (10), i.e., we do the transformation 7Tj — >■ 7Tj + 7Tj and neglect the background 



part of the Lagrangian: 



= - 6/3 d^k^d^j) + 4/3 fc^dpTTidV^d 



TT; 



6^ ] d p 7r i d°7r i d fl dWd x d T 7r j . (11) 
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Requiring the self-accelerating background to be ghost free gives rise to /3 > 0, so the 
conditions for a ghost free self-accelerating solution are 



(5 > and a < and Ek = H > 0. (12) 

Therefore, in the SO(N) (fundamental) multi-galileon theory, when the self-accelerating 
branch is ghost free, the Minkowski branch is inevitably haunted by ghosts, and vice versa. 
Also, we see that there is just one canonical kinetic term on the self-accelerating background, 
while on the Minkowski background there are N of them 1 . 

All of these would become apparent from a point of view of spontaneous symmetry break- 
ing. To facilitate this approach, we would like to utilize the action polynomial introduced 
in [13] 2 : 

f d 4 x L 

L{k) = -4 J, 77 (13) 
K J J d 4 x x^x* v ; 

= ak% + ^(k%) 2 , (14) 

where 7r is evaluated at —kiX^/A. By explicit calculation [13], we have shown that the 
extrema of L(k) correspond to cosmic background solutions; also, the coefficient matrix of 
the canonical kinetic terms of the iV-galileon about a background (ki = is equal to the 
Hessian of L(k) about the background: 



K ij (k) = H v (k)- i)2La) 



dk l dki 



(15) 



k=k 



meaning among the extrema only the (local) minima are ghost free ones. These properties of 
L(k) allow us to treat L(k) as some kind of effective potential in finding ghost free vacua. As 
an aside, note that in canonical field theories such as a scalar field theory, the Hamiltonian 
of the theory provides an energy function to minimize to find stable vacua. However, due 
to their non-trivial vacuum configurations and higher derivative nature, the Hamiltonian 
formulation of multi-galileon theories does not give rise to such a clear energy function for 



the background configuration n = —k^x^x^ /A; see Section IV for details. 

Now, we can easily recover the results of the SO(N) multi-galileon vacuum solutions 
using L(k). The extrema of L(k) give rise to the Minkowski background fcj = and the 
self-accelerating background k^kj = —a/3/3. The background k^kj = —a/3/3 is a minimum 
of L(k) only if a < and /3 > 0. Also, since the continuum k^kj = —a/3/3 is a minimum, 
topologically ki = can not be a minimum, thus for the same set of parameters only one 
of the two backgrounds can be stable. The Hessian of L(k) about the self-accelerating 
background is given by Kijik) = 12f3kikj, which has only one non-zero eigenvalue, so there 
is just one canonical kinetic term on this background. Indeed, we might visualize L(k) with 



1 The same result was also reached in |15j as we were preparing this paper. 

2 Note that here we define a slightly different L(k) from that defined in our previous paper. This is because 
here we write ~ — dirdTrddir..., while in [T3] we use C' n ~ TrddTnrddTT.... These two forms are related 
by integration by parts in the action, so they are physically equivalent. However, when tt is evaluated at 
7r = — kiX^x^ /A, total derivatives also give rise to terms proportional to X jj^X^ • so they differ by a factor 
of —2, i.e., C'^ — —2C n at 7r = —kiX^x^/A. 
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FIG. 1: The Mexican hat shape of the action polynomial L(k) of the SO(N) multi-galileon (plotted 
for the case of 50(2)). The vacuum rests on the trough, spontaneously breaking SO(N) to SO(N— 
1), therefore fluctuations along the N — 1 flat directions of the trough do not have canonical kinetic 
terms. 

a ghost free self-accelerating background by a "Mexican hat" (Fig. [I]). The trough of this 
Mexican hat is an (N — l)-sphere, respecting SO(N). An (iV — l)-sphere (or SO(N)) has 
N(N — 1)/2 independent rotational symmetries. The vacuum solution occupies one point on 
the trough and thus only respects an SO(N— 1) subgroup, which leaves a sub (iV— 2)-sphere 
still rotational symmetric and breaks N — 1 rotational symmetries. Only the radial direction 
around the trough accommodates non-trivial "oscillations", reflecting the presence of only 
one canonical kinetic term. The N — 1 flat directions represent the loss of iV — 1 canonical 
kinetic terms. 



B. General Proof 

This is of course reminiscent of Goldstone's theorem for a canonical scalar field theory 
with a potential. Here we are able to prove an analogous theorem for a symmetric multi- 
galileon theory with an arbitrary internal group that the number of canonical kinetic terms 
that are lost is equal to the number of spontaneously broken symmetries, which in turn 
equals the dimension of the total symmetry group minus that of the unbroken subgroup. 
Again it is sufficient to use the action polynomial L(k) to prove this. 

Let k = k{ be a (local) minimum of L(k), so it is a sensible background to expand the 
theory. Since k = ki is a minimum, Kij(ki) should only have non-negative eigenvalues. 
The eigenvectors of positive eigenvalues correspond to the canonical kinetic terms, while 
the eigenvectors of zero eigenvalues correspond to the degrees of freedom without canonical 
kinetic terms. 

To prove the theorem, we must show that every spontaneously broken symmetry gives 
rise to an independent zero-eigenvalued eigenvector. Under an infinitesimal group action, 
for the configuration 7r' = — A;*x M x M /4, we have 

k'^k' + e A%k), (16) 



where e is an infinitesimal. Since C n is invariant under a group transformation, from (13), 
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we infer that L(k) is also invariant. So we have 



L(k) = L(k + e A(k)) 



m+e^A^k), 



dk l 



which leads to the identity 



dL(k) 
dk l 



A*{k) = 0. 



(17) 



(18) 



Differentiating it with respect to k % and evaluating it at the vacuum of the theory (ki = kj) 



gives 



Kij(k)A j (k) = 0, 



(19) 



where Kij(k) is the coefficient matrix of the canonical kinetic terms, as defined in (15) 



Now, if the transformation (16) belongs to the unbroken subgroup, the vacuum ki 



ki is 
0. If 

the transformation (|16|) belongs to a spontaneously broken symmetry, the vacuum is changed 



invariant under the transformation and the relation (19) is trivial as we have A J (k) 



along the flat directions of the continuous minimum of L(k) and so we have A J (k) 0. In 
this case, Kij(k) has a zero eigenvalue and the eigenvector A^(k), or A J '(7r), is the degree of 
freedom that loses its canonical kinetic term. 



C. Implications 

In multi-galileon theories, due to the presence of higher order kinetic terms, absence of 
a canonical kinetic term does not necessarily mean loss of a dynamical degree of freedom. 
Taking the SO(N) multi-galileon theory for example, by integration by parts the cubic term 



of the Lagrangian above the self-accelerating background (11) can be cast as 



4 3) = - W«9 a (/^,) - %(3d a d a {¥ui) jr% + 4/3 S^d^d^ d c d d (k^,), (20) 

where a,b,c,d are spatial indices (rather than group indices), and the theory has N cubic 
kinetic terms. The conjugate momenta of 7r l (x, t) is non-vanishing and the canonical phase 
space is non-trivial for all the N degrees of freedom. So there are still iV apparent dynamical 
degrees of freedom on the self-accelerating background. 

However, since a mode without a canonical kinetic term can be regarded as infinitely 
strongly coupled and Vainshtein mechanism takes effect in galileon models, these modes 
would self-screen themselves from the others. We shall demonstrate this schematically. 
Suppose 7Ti loses its canonical kinetic term around the self-accelerating vacuum and consider 
a slightly different background where the Lagrangian is given by 

M 2 

C ~ e 2 Mpdn 1 d7r 1 + — £07ri07rid07n + tqT + ... , (21) 
M 2 

with ... standing for other interactions and modes. To see the genuine dynamics of this 
mode, we canonically normalizing it, which gives rise to 

1 M 2 1 
C ~ Mpdiridiri + -—^dTTid^ddTT! + -tx{T + ... . (22) 
e A M l e 
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We can recover the perturbative Lagrangian around the vacuum by taking the limit e — > 0, 
where we can clearly see that Ti\ is infinitely strongly coupled. Now, we can calculate that 
the Vainshtein radius of a spherical source (M s ) for tti (see e.g. [TB"]): 



4 1 

3 



= [jM- (23) 

It goes to infinity when e goes to 0, meaning this mode would be self-screened at infinitely 
large distances, and thus is effectively non-dynamical on the vacuum, at least in terms of 
weak gravitational interactions. Nevertheless, as mentioned above, although some modes in 
the galileon multiplet lose their canonical kinetic terms on a self-accelerating vacuum, these 
modes can re-acquire their quadratic kinetic terms on backgrounds with matter sources. 
Therefore, around a generic background such as in the solar system, these modes are indeed 
not strongly coupled. As an aside, if there is a cosmological constant, the multi-galileon 
internal symmetry will be explicitly broken, in which case there is generally no loss of 
canonical kinetic terms. For the example of 50(2) fundamental representation, the action 
polynomial is deformed to be a tilted Mexican hat where there is only an unique minimum. 

When calculating the leading corrections to GR, thanks to Vainshtein effect, we might 
simply exclude these inert modes. So the spontaneous symmetry breaking and the sub- 
sequent freeze-out of some dynamical modes could be reflected in tests of modification to 
gravity force, as the leading corrections are encoded in the canonical kinetic terms. We 
still take the SO(N) multi-galileon for example. First, note that the one particle exchange 
amplitude between two conserved sources T^ v and T' in GR schematically is given by 

• 4 -~-s?( T "o T -4 T o r )' < 24 > 

For simplicity, we assuming fcj ~ k, so we have Nk 2 ~ —a/3/3. When the vacuum is 



spontaneously broken and rests on the self- accelerating branch (10), from (11) we can see 
that the SO(N) multi-galileon gives rise to a leading correction 

N 1 

5 Asa = A-A gr ~ ^T-T'. (25) 

This is compared to the case without spontaneous symmetry breaking, when the leading 
correction on the Minkowski branch ^ is given by 

SA M ~ --T^T'. (26) 

On the other hand, when testing the multi-galileon modification to gravitational force 
upto leading order, we have to deal with observational degeneracy between multi-galileon 
theories with different internal symmetries and different choices of vacuum branches. Again 
taking the SO(N) example and assuming ki ~ k, the leading correction from the SO(N) 



multi-galileon on the self-accelerating branch (25) is the same as that from a multi-galileon 



theory without internal symmetries and with canonical kinetic terms — a(dirid7ri + ... + 
dTT N / 2 dn N / 2 ), provided N is an even number. 
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IV. HAMILTONIAN FORMULATION OF MULTI-GALILEON THEORIES 



In this section, we deviate from our main plot of the paper and briefly introduce a subplot: 
the Hamiltonian approach of multi-galileon theories. First, we derive the Hamiltonian for a 
general multi-galileon theory with or without internal symmetries. 

As the Lagrangian of a multi-galileon theory contains terms with more than 2 spacetime 
derivatives, one might expect the Hamiltonian formulation of a multi-galileon theory should 
involve Ostrogradski's prescription for high order derivative theories (see for example [21]). 
However, a bell should be certainly rung to this naive thinking once we notice the fact that 
the equations of motion of a multi-galileon theory has only second order derivatives. We will 
see that a general multi-galileon Lagrangian can be cast to have only up to first order time 
derivatives. A general multi-galileon theory without a tadpole term can be written as [H] 



a' 



d, 2 ir n d^ l2 d^H l3 ...d, n d^ ln 



(27) 



where S^'"^ 71 -, = (n — l)\8fi£...6£\, ii,...,i n label different galileons (not necessarily inter- 
nal group indices) and summation over repeated ik is understood. a n "- ln are free pa- 
rameters of the theory, and can be chosen as symmetric in exchanging the indices since 

^\yz''vn\ ^2 ff n^ 27r «2^3^ 37r «3---^n^ I/n7r jn can De made symmetric in exchanging the galileon 
indices by integration by parts. 

To see what the derivative structure is, we should unfold the anti-symmetrisation. Since 
the Hamiltonian formulation only requires the knowledge of the time derivative structure, 
we only need to separate the time derivatives from the spatial ones. A useful relation for 
the separation is 



r fJ2---fJ n rpV2..M n 
[u 2 ...U n ] fJ.2---fJ.-n. 



C' 

[02 ■■ 



rpb2-..b„ _i_ \ A \ Sj 0-2---0-TL rpbi—bn 
A a.2...a„ ' / j / j u [b 2 ..-b n ] «2-in 
i=2 3=2 



►tl 
>1 2 



(28) 



where is an arbitrary tensor, t\ and t2 are time indices, and and hi are spatial 



indices. The double summation is over replacement of one up spatial index with ti and one 
down spatial index with t 2 , so there are (n — l) 2 terms with time derivatives. Applying this 



formula to (]27|) and repeatedly integrating by parts, we can see that for n-th order a term 

7r il Tr i2 d a p b3 TT i3 ...d a ,d K 



with <$■;;££ 



gives rise to a 



, ca 3 ...o 

°lb 3 -b- 



'7T,' 



while a term with S"^ 1 ^' 



only gives rise to half of that, with all the other terms cancelling each other. Therefore the 



Lagrangian (27) can be cast as 



n=2 



d a p b ^ la ...d an d K v ln 



(29) 



where C 2 = n(n — l)/2. The appearance of the combinatorial number C\ is what one might 
expect, since the indices i\, ....in are symmetric and so there are C\ ways to pick out two 
7TjS with first order time derivatives. Due to the first order structure in time derivatives, we 
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can simply take 7Ti(x, t) canonical coordinates and define the conjugate momenta as 



l (x,t) 



07Ti(x,t) 
5 

2 £ n'- * l2 d a p b ^...d an d b ^ n . 



n=2 



Defining the matrix 



n=2 



we can reverse ( 30 ) and get 



7T, = (M~ X )ij$ 



(30) 



(31) 



(32) 



To get the Hamiltonian of the multi-galileon theory, we perform the Legendre transformation 



H w = d 3 x 



TliCj) 1 - C w 



d 3 x^, 



(33) 



where the Hamiltonian density is given by 



n=2 
1 



(34) 



n=2 



Now, we would like to know what the Hamiltonian looks like for the vacuum configuration 
7Ti = -kiX^/A: 



H n = ] [ d 3 x \(t 2 + x 2 )L^(A;) + (3t 2 + 2 -x 2 )L^\k) + (6t 2 + \x 2 )L^(k) 
4 J 3 3 



+ (10t 2 )L (5) (*0 



(35) 



where L^'(k) are the z-th order terms of the action polynomial L{k). In a canonical field 
theory with a constant field background, since the Hamiltonian is an (infinitely) extensive 
quantity, we can divide the Hamiltonian by the volume of the spacetime to extract an energy 
function of the constant field, which can be minimized to find the vacua of the theory. Here 
we find the same procedure is not applicable to a multi-galileon modified gravity theory, as 



we can see from (35) that the "volume factor" is different for different orders of fcj. This 



of course originates from the high derivative nature of multi-galileon theories and the non- 
trivial background configuration 7Tj = —kiX^x^/A. Note that for a multi-galileon action with 
the configuration 7Tj = —kiX^x^/A, a total derivative (say, dt(iiid t 7i2did , "7i 3 )) will actually give 
rise to nontrivial contribution (— 3/ci/c2^3(3t 2 — x 2 )/16). Indeed, from (27) to (29) we have 



performed a series of integration by parts and neglected the subsequent total derivatives, 



which is responsible for the different "volume factors" in (35) 
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We also note that the Hamiltonian density (34) (hence the Hamiltonian) is generally 
unbounded from below, i.e., the Hamiltonian density can be arbitrarily lowered by choosing 
suitable initial field configurations. This is due to the presence of higher than quadratic order 
multi-galileon terms and because there are terms where the first derivatives of galileon fields 
are not in "squared" forms (e.g., ^i7ri(<9 a <9 a 7r 2 ) 2 is "squared", but d a ii 1 d a dbTCid b Tiid c d c Tii and 
7ri7r 2 <9 a <9 a 7T3 are not.). Since we know the galileon models define a conventional Cauchy prob- 
lem, the galileon fields and their first derivatives can be arbitrarily chosen. By making the 
first derivatives of galileons increasingly steep, we can lower the Hamiltonian density arbi- 
trarily. Note that, even for the background configuration 7Tj = —kiX^x^/A, the Hamiltonian 
density at a fixed spacetime point is not bounded below if the highest order of galileon terms 
is odd. The perturbative Hamiltonian above some self-accelerating background (fcj = kj) 
can also be cast in the form (33) with the parameters a 11 '"'" replaced by a new set of pa- 
rameters /3 n - in (^) (as polynomials of ki) (see e.g. [13]), so it is also unbounded below. In a 
fundamental theory, this of course signals instabilities. However, the multi-galileon modified 
gravity is only supposed to be the decoupling limit of some underlying full theory, so one 
should really check whether the Hamiltonian of the underlying full theory is well behaved 
or not. The underlying theory presumably has 4D diffeomorphism invariance, so the cor- 
responding naive 4D Hamiltonian (excluding the part from extra dimensions) is tuned to 
zero by 4 constraint equations, similar to that in GR. A useful 4D Hamiltonian arises when 
the theory is "deparameterized" [2S], but from the experiences in GR, even checking the 
positivity of the background solution could be nontrivial 3 . 

On the other hand, due to the derivative structure of the multi-galileon theories, the 
most negative Hamiltonian value is achieved by setting the gradients close to the cutoff of 
the theory, i.e., d ~ A cuto fr- This kind of being unbounded below pushes the limit of a 
classical theory, as it relies on a small region of the canonical phase space, so one might 
also doubt whether quantum corrections can alther the picture. A famous example of this is 
the Hydrogen atom: The classical Coulomb potential for this system (— e 2 /r) can be made 
arbitrarily negative by placing the electron close to the nucleus, but the Hydrogen atom is 
stable upon quantisation of electrodynamics. 



V. CONCLUSION 



We have coupled the multi-galileon theory with internal symmetries studied in [T4] to 
conventional General Relativity (GR) and proposed it as a modified gravity theory in the 
decoupling limit where the multi-galileon modifies GR only by mixing with the transverse 
graviton. We have discussed the phenomenon of spontaneous symmetry breaking of these 
theories on (classical) self-accelerating vacua. We point out that, similar to that in canonical 
scalar field theories, the pattern of the symmetry breaking is governed by a new version of 
Goldstone's theorem that for every broken continuous symmetry the theory loses a canonical 
kinetic term. Note that as the energy-nomentum tensor T^ u by definition vanishes in the 
self-accelerating vacuum, this theorem is largely insensitive to the coupling to GR. But we do 
assume the background configuration of the multi-galileon is given by 7Tj = —kiX^x^/i. We 



Nevertheless, for the case where the multi-galileon Hamiltonian density ( 33 ) is included in the constraint 
equations, the multi-galileon Hamiltonian (33) being unbounded below is irrelevant to the stability issue 
of the full theory. 
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have also discussed implications of this theorem. In particular, we suggest that the mode 
that loses its canonical kinetic term, although apparently non-trivial in the phase space, 
becomes inert due to Vainshtein mechanism. This would lead to different modification to 
gravitational force, compared to what one would naively expect from the Lagragian with 
the broken vacuum hidden. Also, there would be degeneracy among multi-galileon theories 
with different internal symmetries and different choices of vacuum branches. 

We have also derived the Hamiltonian of a general multi-galileon theory. We find the 
Hamiltonian with the configuration 7Tj = — kiX^x^/A does not give rise to a good "effective 
potential" to minimize to find the background solution. Besides, we find the Hamiltonian 
is not bounded below because of the higher order multi-galileon terms. We speculate this 
pathology might arise from the decoupling limit or the classical nature of multi-galileon 
theories and argue that the underlying full theory for the multi-galileon or even its quantum 
version should be investigated to decide whether this is a real problem or not. There are a 
few attempts to put the galileon description in a more formal framework [5j [TJ 19] [15] , and 
it is interesting to see whether the vacuum Hamiltonian in these models is bounded below, 
which we leave for future work. 
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